*Figure A6
xtset
preserve
tempname myresults
local i=0
postfile `myresults' var beta se using myresults.dta, replace
foreach v of varlist sum_shame sum_resource sum_coord {
local dv agg_implem_wo 
local iv lag_runn lag_l* infant_wb lag_constraint merged_gdp 
local i=`i'+1
	foreach model in ///
	"xtreg `dv' `v' `iv', fe cluster(pam_caseid)" ///
	"reg `dv' `v' `iv', cluster(pam_caseid)" ///
	"xtreg `dv' `v' `iv', cluster(pam_caseid)" ///
	"reg D.`dv' D.`v' D.(`iv'), cluster(pam_caseid)" ///
	"xtreg `dv' `v', fe cluster(pam_caseid) "{
		`model'
		matrix b=e(b)
		matrix se=e(V)
		scalar beta=b[1,1]
		scalar ses=se[1,1]
		post `myresults' (`i') (beta) (ses)
		
		}
	
	
}
postclose `myresults'
restore

use "myresults.dta", clear
replace se=sqrt(se)
g ub=beta+1.96*se
g lb=beta-1.96*se
egen model=fill(1 2 3 4 5 1 2 3 4 5)
label define model1 1 "Fixed Effects" 2 "OLS" 3 "Random Effects" 4 "First Differenced" 5 "Bivariate FE"
label values model model1
label define varl 1 "Shaming" 2 "Funding" 3 "Coordinating"
label values var varl

foreach v in 1 2 3 {
twoway rcap ub lb model if var==`v', lcolor(black) || (scatter beta model, mcolor(black)) if var==`v', legend(off) ///
xtitle("Specification") xlab (1(1)4) xlabel(1 "Fixed Effects" 2 "OLS" 3 "Random Effects" 4 "First Differenced" 5 "Bivariate FE", angle(45)) ///
ytitle(Effect of `v') yscale(r(0 14)) ylabel(0(2)14) graphregion(color(white)) saving(`v'.gph, replace)
}

graph combine 1.gph 2.gph 3.gph, graphregion(color(white))
graph export "Figure A6.pdf", replace

use "Remaining Seized Replication Data.dta", clear
